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Abstract. Many past attempts to kill MOND have only strengthened the theory. Better data 
on galaxy velocity curves clearly favor MOND (without fine-tuning) over cold dark matter. 
The usual critism on the incompleteness of classical MOND has spurred a Modified Relativity 
(MR) by Bekenstein. After outlining cosmology and lensing in MOND, we review MOND on 
small scales. We point out some potential problems of MOND in two-body relaxation and 
tidal truncation. We argue that the tidal field in any MOND-like gravity theory predicts that 
the Roche lobe sizes of a binary system are simply proportional to the binary baryonic mass 
ratio to the power 1/3. An immediate application of this result is that the tidal field and tidal 
truncation radii of million-star globular clusters and million-star dwarf galaxies (e.g., the Milky 
Way satellites NGC2419 and Carina) would be very similar because of the one-to-one relation 
between gravity and baryon distribution. This prediction appears, however, inconsistent with 
the fact that all globulars are truncated to much smaller sizes than all dwarf galaxies. Whether 
tide is uniquely determined by baryons can also be used to falsify any MOND-like gravity theory, 
whether classical or relativistic. 
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1. Introduction 

A large amount of astronomical data suggests that the GM/r 2 Newtonian gravity 
from the baryonic material in galaxies fails to explain the large accelerations implied by 
motions in spiral and elliptical galaxies. However, Zwicky's proposal to introduce unseen 
gravitating (dark) mass with a very flexible density distribution wherever needed for the 
motions of galaxies and stars has been troubled by a lack of experimental confirmation of 
its elementary particle counterparts for more than 60 years. Evidence is also summarized 
in recent articles about the difficulties facing dark matter (DM) theories on galaxy scales 
(Spergel & Steinhardt 2000, Ostriker & Steinhardt 2003). Surprisingly, an alternative 
gravity theory with Modified Newtonian Dynamics (MOND, Milgrom 1983) has been 
doing very well exactly where the dark matter theory is doing poorly. The predictive 
power of this 20-year-old classical theory with virtually no free parameters (Bekenstein 
& Milgrom 1984) is recently highlighted by the astonishingly good fits to contemporary 
kinematic data of a wide variety of high and low surface brightness spiral and elliptical 
galaxies; even the fine details of the ups and downs of velocity curves are elegently repro- 
duced without fine tuning of the baryonic model (Sanders & McGaugh 2002, Milgrom 
& Sanders 2003). Originally it was proposed that galaxy rotation curves could be fit by 
motions in a modified gravity g, which is significantly stronger than Newtonian gravity 
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in the weak regime defined by a gravitational field energy density 

where G is the gravitational constant, and pb{0)c 2 is the present universal mean baryon 
energy density (commonly the weak regime is defined as where gravity g ^ ao ~ 1.2 x 
10~ 8 cmsec~ 2 ). Far away from the baryons, the modified gravity g satisfies 

En = MA 



G(g) G - ^ \Rtf 

where gN is the vector sum of Newtonian gravity (distance-inverse-squared attractive 
force) of all baryonic particles with masses Mj at distances Ri, where the effective gravi- 
tational constant G can be chosen as, e.g., -^L = 1 — exp — ( — ) , where a = 1 gives 

G(g) l \ a o J J 

the MOND gravity g. We get a MOND-like theory for any positive a. In these theories, 
the only matter that matters is the luminous (baryonic) matter, and the gravitational 
field is a unique function of the luminous matter distribution. 

These alternative theories based on modifying the law of gravity in a purely baryonic 
universe have been gaining ground rapidly. Milgrom's empirical formula is a simplified 
treatment of a classical theory with a curl-free gravitational field by Bekenstein & Mil- 
grom (1984), which is now extendable to a Modified Relativistic (MR) theory (Bekenstein 
2004) which passes standard and cosmological tests used to check General Relativity 
(GR). This recent breakthrough by Bekenstein transforms the empirical MOND formula 
to a falsifiable theory with respectable theoretical underpinnings comparable almost to 
the rival dark matter theory. Here we examine various aspects of the MOND theory after 
outlining cosmology and lensing in the relativistic MOND. 

1.1. Cosmology and Lensing 

Bekenstein (2004) has shown that MONDian cosmology and lensing follow the standard 
GR formulation. The physical metric of the FRW universe in MOND is written as, 

9flv dx^dx u = c 2 dt 2 - R{t) 2 [d X 2 + f(x? (d0 2 + sin 2 8d<j> 2 )] (1.3) 

where f(x) = sinx,x,sinh% for closed, flat and open universe as in FRW GR, and the 
physical scale factor is R(t) = Ro/(l + z). Note that we have removed the tilde sign that 
Bekenstein uses to denote physical coordinates, if we introduce the vaccum energy and 
neglect the scale field <j> as in Bekenstein, then the MONDian universe expands with a 
Hubble parameter H{z(t)) = d ln ^^ given approximately by 

^fi- « n m (i + zf + n,.(i + zf + q k (i + z) 2 + n A ,n K = i- n m -si r - n A) (i.4) 

where fl r ~ 10~ 4 /i~ 2 is normalised by the 2.7K CMB, and Q m ~ 0.05h~ 2 from all 
observable baryons, and Q,k is the curvature term. The matter-radiation equality happens 
at redshift of about f2 m /f2 r ~ 1000; this is much lower than in CDM models. 

To do lensing, Bekenstein shows (in eq. 116) that for small perturbations of the metric 
by a potential <&, light travels on geodesic given by 

2<T> 2<J> 
ds 2 = g^dx^dx" w -(1 + -^)c 2 dt 2 + (1 - -j)dl 2 = 0, (1.5) 
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where all is the proper length. Hence the light travel time 
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same as in Einstein universe. The gravititional potential satisfies approximately the 
MOND eq. 
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where p is the gravitational mass density of baryonic matter, and the fi is the usual 
dimensionless function of the rescaled gravity introduced by Milgrom. 

The deflection for a light array of the closest approach (roughly impact parmameter, 
physical length) b is calculated by 

2gj_olt 



a 



(1.8) 



where g± = g{r) ■ £ is the perpendicular component of the gravity g(r) at radius r = 
\J (ct) 2 + b 2 , and ct is the path length along the line of sight counted from the point of 
closest approach. Hence the deflection angle 

2g(r)alt b 



J — ( 



C y/(ct) 2 + b 2 

This is essentially the same as eq. (109) of Bekenstein. 



(1.9) 



2. Two-body Relaxation Time in MOND 

As pointed by Ciotti & Binncy (2004) , relaxation in MOND is a non-trivial calculation 
because it is a long-range gravity with a divergent potential at large distance. However, by 
and large, MOND gravity is effectively a stronger-than-Newtonian gravity. This suggests 
that we can rescale Newtonian gravity by a factor ^. To balance this stronger gravity, 
a bound object of certain size r acquires a higher internal velocity (than Newtonian 
case with pure baryons), hence shorter dynamical time, and faster relaxation since the 
relaxation time is proportional to dynamical time with 

Nr 

t x OC Nt dyn ~ (2.1) 

where C ~ 71n./V ~ 50 is a large constant related to the Coulomb logarithm, and the 
dynamical time 

t*yn=(y) -lMyr-— — (2.2) 

where r is the half-mass radius, and V is characteristic circular velocity. So for TV <~ 10 4 ~ 5 
globular cluster, it relaxs on time scales of 100-1000 dynamical times. Typical dynamical 
time is of order Myr, hence a globular is generally relaxed in a fraction of the Hubble 
time. In MOND the circular velocity is generally given by 

/ TUT \ 1/4 

V = (GMa f 4 = S.Skms- 1 (^i^J (2-3) 

if in dcep-MOND where r ^ = 3.5 pc^Jj^^. 

Relaxation drives a globular to expand. Star clusters might have been born very com- 
pact, and in the strong-gravity regime. However, fast relaxation means that globular 
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Figure 1. shows the relaxation-driven expansion of the half- mass radius of an isolated globular 
cluster of N = 2 x 10 4 ,6 x 10 4 , 1 x 10 5 (curves from top to bottom) half-solar- mass stars with 
a pre-core-collapse size of lpc (marked by the blue circles). The change from strong to weak 
gravity is indicated by the change from thin lines to thick lines. Observed globulars typically 
occupy the green shaded zone in the age vs. size plane. 



clusters cannot stay small, will eventually expand and enter the weak regime. The relax- 
ation time can be expressed in terms of some kind of rate of expansion per relaxation 
time 



where we have assumed C ~ 71niV ~ 50, and have used the circular velocity in the 
MOND regime. The above implies that a cluster should expand to a size of ~ 175pc over 
a Hubble time in MOND. 

To illustrate this, we give a few examples in Fig.^ Consider, for example, a hypothet- 
ical cluster born fully isolated with N — 2 x 10 4 half-solar-mass stars inside a half-mass 
radius of lpc. The cluster is initially dense and in the strong gravity regime. Rapid re- 
laxation leads to core collapse after about to = 0.4 Gyrs where we take to as 10 times the 
initial relaxation time t x - Afterwards the cluster expands as t 2 ^ 3 in strong gravity till rh 
reaches about 3pc at about t = 1.4 Gyr, after which the cluster enters the weak gravity 
regime, Vh expands linearly with time according to MOND, and reaches a size = 45pc 
after a Hubble time. Systems with larger N grows slower (Fig. QJ. Nevertheless, clusters 
with N = 2 x 10 4 to 1 x 10 5 half-solar-mass stars all grow to a half-mass size 12 — 45 pc, 
too big to be consistent with observed globular clusters. 




(2.4) 
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Figure 2. shows the rescaled Roche lobes (contours of the effective potential) in the equatorial 
plane of a hypothetical isolated Earth-Sun binary mass ratio 3 x 10 -6 in the strong gravity 
regime (say with the separation 1AU, thin blue lines), and in the weak gravity regime (say with 
separation lpc, and GMq/1 pc 2 ~ O.OOlao <C flo> shaded areas). The Earth is at origin and the 
Sun is at rescaled unit length to the right. The inner Lagrangian point is a saddle point between 
the Earth and the Sun, which is slightly further away from the Earth in the deep-MOND regime 
than in strong gravity regime. 



3. Tidal truncations or Roche lobes of binary systems 

First we note that at large distances from the baryonic particles Gauss's theorem 
predicts that a MOND-like gravity g with an effective gravitational constant G(g) is 
approximately parallel to the direction of the Newtonian gravity gjv with the amplitude 
satisfying eq. (JIJ, so the gravity g is determined by the total mass of the baryons J2i Mi 
and the mean distance r, independent of the spatial distribution of the baryons. This 
is true even in the rigorous classical counterpart (Bekenstein & Milgrom 1984) of the 
Bekenstein theory. So, e.g., a dense globular cluster and a fluffy dwarf galaxy of 10 5 L Q 
should have comparable baryonic mass, hence the similar gravity field on scales of kpc, 
nearly independent of the internal density profile. Indeed MOND nicely accounts for the 
very small scatter of dynamical properties of structures of similar baryonic distribution 
on galaxy scale (Sanders McGaugh 2002, McGaugh 2000, 2004). However, this "nice 
feature" of baryonic gravity also means that the dynamics of wide binary stars would be 
a scaled-down version of the Antenna-like merging galaxies. 

Theoretically in a binary system, the Roche radius (or tidal radius or virial radius 
when referring to galaxies sometimes) is where the average density of the mass-losing 
satellite is comparable to the average density of the mass-receiving host at some point R 
along the orbit of satellite, i.e., 

Gm GM 

fc 4^7^> I 3 - 1 ) 



with the prefactor k ranging from 1 to 3 in different definitions in the literature (Binney 
& Tremaine 1987); e.g., k = 2 according to von Hoerner (1957). This well-known tidal 
criterion in Newtonian gravity is in fact a result of basic dynamics, independent of the 
force law parameter G(g) and applies to any gravity. It is shown elsewhere (Zhao et al. 
2005) that the prefactor k = 2 rigourously in the deep-MOND regime. The shape of the 
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Figure 3. shows the satellite distance vs. luminosity scatter diagram for the dwarf galaxies 
(labeled with larger font) and globular clusters (labeled with smaller font) orbiting in Milky 
Way's very outer halo (with orbital periods about one quarter of a Hubble time). The sizes of 
the symbols are scaled to the physical sizes of the satellites (as a reference, the faint Eridanius 
cluster at the bottom is lOOpc, and Sextans dwarf is 3000pc). The sizes of all globulars of the 
Milky Way are all drawn as smaller rings inside NGC2419, which is the largest, and most distant 
bright globular cluster. Predicted sizes of the faint Pal3 cluster and Carina dwarf are also shown 
as squares based on their proper motion inferred pericentres in a 10 n Z/Q baryonic Milky Way. 



Roche lobes is somewhat more squashed than the Newtonian counterpart. An example 
of Roche lobes is given in Fig. [3J 

3.1. Application to satellite globular clusters and dwarf galaxies 

The Roche criteria (eq. I3.1[l predicts that objects should have same sizes if they have 
similar baryonic content in similar enviornment. This, however, is inconsistent with basic 
observation data, and we see a large scatter for the truncation radii of outer globular 
and a dozen dwarf galaxy satellites accelerating in the tide and very weak gravity of a 
luminous host galaxy, like the Milky Way (see Fig. [3J • So the observed limiting radii of 
these satellites (Harris 1996, Mateo 1998) are no longer simply and uniquely scaled with 
the binary baryonic mass ratio as predicted by a MOND-like theory (eq. 13.1(1 . 

For example, NGC 2419, the largest cluster of the Milky Way, is as distant as and 
1-2 times as luminous as the (2.5 — 4) x 10 5 L Q dwarfs Carina/Sextans, which has an 
equilibrium of stars extending to 1000-3000pc radii at a distance of about 90kpc from 
the Milky Way. A baryonic universe predicts that NGC 2419 should resist the tide of 
the Milky Way as well as the Sextans dwarf does on 3000pc scale; Sextans shows no 
sign of being disrupted (Grebel et al. 2000). Nevertheless, even this largest globular is 
truncated on 200pc scale. Such large discrepancy does not go away if we take into account 
of different orbits or mass of the satellites, e.g., a baryonic universe would predict that 
the cluster Pal3 (faint 10 4 i© and presently near pericentre at 90 kpc, Dinescu et al. 
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1999) is bigger than the Carina dwarf (presently near apocentre, Piatek et al. 2003), 
but in reality as shown in Fig. |3| Pal3 (Sohn et al. 2003) is much smaller than Carina 
(Walcher et al. 2003). Note that orbits in MOND is not much different from orbits in 
isothermal halo potentials (Read & Moore 2005). 

This example points out a generic problem general to galaxies with resolved satellite 
globulars vs. dwarfs in any baryonic universe including Bekenstein's cosmology: why there 
is a clear size gap between 250pc-400pc and nature forms neither very big globulars nor 
very small dwarfs despite a good mix of the baryonic mass and enviornment variables (or- 
bital phase and pericentric tide) of these satellites. In a baryonic universe "globulars" and 
" dwarf ellipticals" are practically similar assembly of stars described by the King models: 
they are pressure-supported low-dispersion gas-poor low-luminosity elliptical equilibria 
immersed in and sculptured by the tidal field of a common host galaxy. 



Could it be that the globulars were born small and stayed small? It is plausible that 
globulars might have formed from denser and cooler gas clouds than stars in some dwarfs. 
However, a very dense core of a globular is theoretically unstable to stellar encounters and 
must evolve rapidly, and on a Hubble time typically 10-100 percent of the initial mass of 
a star cluster are unbound due to relaxation-driven stellar evaporation (Johnstone 1993). 
Relaxation also drives the half-light radius of a globular to expand as a power law of the 
time t 2 / 3 if we assume globulars are born infinitely compact and dense so we can neglect 
the delay for the core to collapse. Using eq. (1) and Fig. 1 of the numerical simulations 
of (Giersz & Heggie 1996), we estimate the radius containing 90 percent of the mass of 



Relaxation is even faster in MOND (Ciotti & Binney 2004) because of the stronger 
gravity, higher velocity dispersion and faster dynamical time (Baumgardt et al. 2005). 
A detailed calculation of the expansion rate in MOND is given elsewhere (Zhao et al. 
2005). Essentially a globular cluster born very compact in MOND cannot stay compact. 

There are five bright old globular clusters orbiting around the Fornax dwarf elliptical, 
140kpc away from the Milky Way centre. These globulars have a typical baryonic mass 
of lO 4_5 Af enclosed in a limiting radius of 30 — 50 pc with a half-mass radius of 4 — 12 pc 
(using data compiled in Mackey & Gilmore 2003 and Rodgers & Roberts 1994). One of 
clusters, No. 1, has a distorted profile and appears to be tidally disturbed. The clusters 
No. 3,4,5 appear to have excess stars near their limiting radii. The clusters No. 1,2, 5 
are on the outskirts of Fornax, between l-1.5kpc from the Fornax centre. Our model of 
relaxation-driven expansion predicts (cf. Fig. ^) that Fornax No.l (1 x 10 4 L & ) should 
have grown to a half-mass size of about 40pc, and more massive clusters No. 2 and No. 
5 (5 x lO 4 ^©) should grow to 12pc. The predicted sizes are significantly greater than the 
observed half-light radius of 12pc (No.l) and 6pc (No. 2 and No. 5). 

Our model also predicts (cf. eq. 13. 1J) that the Roche lobe of these Fornax clusters 
should have a radius of (650 — 1000) pc if they are on a 140kpc orbit around the Milky 
Way (5 x 1O 1O L0), or a Roche lobe radius of (100 — 130) pc if they are on a 1.5 kpc 
orbit around Fornax (2 x 10 7 L Q ). Both of these predictions are much greater than the 
observed limiting radius of (40-50)pc for all Fornax clusters. 



The limiting size to which a globular can maintain internal equilibrium are most likely 
limitted by the tidal force of the Milky Way in general. This is supported by the obser- 
vation of a characteristic sudden change of outer density profile of most globulars and 



3.2. Application: the sizes of the globular clusters in Fornax 



a typical globular cluster of 10 5 half-solar mass stars expands as lOOpc 




3.3. Discussion 
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Figure 4. shows the central line of sight velocity dispersion of Galactic globular clusters (in 
histogram), and the velocity dispersion of Galactic dwarf galaxies (circles). Data taken from 
Gnedin et al. (2002) and Mateo (1998). 



transient non-spherical outer features, as illustrated by the beautiful stream of Pal 5 found 
by Sloan Digital Sky Survey, and many distant clusters, e.g., NGC5694 and NGC1904, 
when deep observations have been taken along lines of sight with least background con- 
fusion (Leon et al. 2000, Odenkirchen et al. 2003, Grillmair et al. 1995,1996, Lehmann 
& Scholz 1997). These features are most likely due to stars leaking from a boundary set 
by the Roche lobe much like the donor star in a binary evolving off the main sequence 
and transfering mass to its companion. The MOND gravity has a ln(r) infinite potential 
well, and so if in isolation a globular will grow indefinitely due to two-body relaxation. 
So both two-body relaxation and violent relaxation tend to smooth out any sharp initial 
boundaries. It is hard to imagine a non-tidal mechanism to sustain sudden changes of 
the density profiles. A globular must have a finite potential well or an outer truncation 
also because large amount of recycled gas from aging and evolved stars must flow out 
of the tidal radii of a globular as old globulars keep less than 1M gas inside (Spergel 
1991). 

The cores of many Galactic globulars are as hot as dwarf galaxies with central velocity 
dispersion ^ 8km/s (cf. Fig. and some dwarf galaxies, e.g., Sextans (Kleyna et al. 
2004) are colder than most globulars. An isolated object enters weak regime if bigger than 
■J GM/ao ~ 3.5 pcy/ Af/10 4 Af Q , hence dwarf galaxies, and outer globulars like NGC2419 
and Palomar clusters are all partly or entirely in the weak gravity regime. Nature some- 
times put star clusters inside a small galaxy (as in Fornax, WLM and Sagittarius dwarf 
galaxies) , but never vice versa. In fact there are evidences that some globulars are stripped 
out from the envelopes or carved from the cores of dwarf galaxies by galactic tide, which 
explains the similar orbits and smooth transition of baryonic mass of these two popula- 
tions of satellites. It seems that while MOND-like gravity theories can account rotation 
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curves of galaxy-scale structures without fine tuning, they do not genericly come with 
an important ingredient or sub-galaxy scale parameter that distinguishes dwarfs from 
globulars. Having dark matter in dwarves but not in globulars is perhaps one way out of 
the dilemma of a baryonic universe. 
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